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Abstract 

Previous findings of Middle East Respiratory Syndrome coronavirus (MERS-CoV)-related viruses in bats, and the ability 
of Tylonycteris- BatCoV HKU4 spike protein to utilize MERS-CoV receptor, human dipeptidyl peptidase 4 hDPP4, suggest 
a bat ancestral origin of MERS-CoV. We developed 12 primary bat cell lines from seven bat species, including 
Tylonycteris pachypus , Pipistrellus abramus and Rhinolophus sinicus (hosts of Ty/onycfer/s-BatCoV EHKU4, Pipistrellus- 
BatCoV HKU5, and SARS-related-CoV respectively), and tested their susceptibilities to MERS-CoVs, SARS-CoV, and 
human coronavirus 229E (HCoV-229E). Five cell lines, including P. abramus and R. sinicus but not T. pachypus cells, were 
susceptible to human MERS-CoV EMC/2012. However, three tested camel MERS-CoV strains showed different 
infectivities, with only two strains capable of infecting three and one cell lines respectively. SARS-CoV can only 
replicate in R. sinicus cells, while HCoV-229E cannot replicate in any bat cells. Bat dipeptidyl peptidase 4 (DPP4) 
sequences were closely related to those of human and non-human primates but distinct from dromedary 
DPP4 sequence. Critical residues for binding to MERS-CoV spike protein were mostly conserved in bat DPP4. DPP4 was 
expressed in the five bat cells susceptible to MERS-CoV, with significantly higher mRNA expression levels than those in 
non-susceptible cells {P = 0.0174), supporting that DPP4 expression is critical for MERS-CoV infection in bats. However, 
overexpression of T. pachypus DPP4 failed to confer MERS-CoV susceptibility in T. pachypus cells, suggesting other 
cellular factors in determining viral replication. The broad cellular tropism of MERS-CoV should prompt further 
exploration of host diversity of related viruses to identify its ancestral origin. 


Introduction 

Coronaviruses (CoVs) are important pathogens in ani¬ 
mals and humans, responsible for a variety of respiratory, 
enteric, hepatic, and neurological diseases. They are now 
classified into four genera, Alphacoronavirus , 
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Betacoronavirus , Gammacoronavirus , and Deltacor- 
onavirus, with Betacoronavirus further divided into 
lineages A to D _4 . Humans are infected by six CoVs, 
including human CoV 229E (HCoV-229E) and human 
CoVNL63 (HCoV-NL63) belonging to Alphacoronavirus; 
human CoV OC43 (HCoV-OC43) and human CoV 
HKU1 (HCoV HKU1) belonging to Betacoronavirus 
lineage A; Severe Acute Respiratory Syndrome-related CoV 
(SARSr-CoV) belonging to Betacoronavirus lineage B; and 
Middle East Respiratory Syndrome CoV (MERS-CoV) 
belonging to Betacoronavirus lineage C' ~ 12 . The emer¬ 
gence potential of CoVs is believed to be related to their 
tendency for mutation and recombination, allowing the 
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generation of new viruses being able to adapt to new 

hosts 3 ' 13 - 18 . 

Bats are an important reservoir of alphacoronaviruses 
and betacoronaviruses, which may jump to other animals 
or humans to cause new epidemics ’ . For example, 
SARS-CoV is likely a recombinant virus originated from 
horseshoe bats as the primary reservoir and palm civet as 
the intermediate host 6,20_25 . Since the SARS epidemic, 
numerous other novel Co Vs from humans or animals 
have been discovered ’ , allowing a better under¬ 

standing of the evolutionary origin of emerging CoVs. 

Although dromedary camels are now known to be the 
immediate animal source of the recent MERS epidemic, 
the evolutionary origin of MERS-Co V remains 
obscure ’ . When the virus was first discovered, it was 
found to be closely related to Tylonycteris bat CoV HKU4 
(Ty-BatCoV HKU4) and Pipistrellus bat CoV HKU5 (Pi- 
BatCoV HKU5) previously discovered in lesser bamboo bat 
(Tylonycteris p achy pus) and Japanese pipistrelle (. Tipis - 

Q 1 1 QQ Q A 

trellus abramus) respectively in Hong Kong ’ ’ ’ . Five 
other lineage C betacoronaviruses closely related to MERS- 
CoV were subsequently detected in bats, including BtVs- 
BetaCoV/SC2013, Hypsugo BatCoV HKU25 from China 
bats, and coronavirus Neoromicia/PML-PHEl/RSA/2011 
(NeoCoV), BtCoVNeo5038/I<ZN/RSA/2015, and BatCoV 
PREDICT/PDF-2180 from African bats _o9 . Besides bats, 
a lineage C betacoronavirus, Erinaceus CoV VMC/DEU, 
subsequently defined as a novel species, Hedgehog cor¬ 
onavirus 1, has also been discovered in European hedge¬ 
hogs, a group of animals being phylogenetically closely 
related to bats °. While none of these animal viruses 
represent the immediate ancestor of MERS-CoV, the spike 
protein of Ty-BatCoV HKU4 was most closely related to 
that of MERS-CoV, and was shown to utilize the MERS- 
CoV receptor, human dipeptidyl peptidase 4 (hDPP4) or 
CD26, for cell entry 1,42 Previous studies have also shown 
that MERS-CoV was able to infect bat cell lines and 
Jamaican fruit bats 3 ’ 44 . These findings suggest that bats 
may be the primary host of the ancestor of MERS-CoV. 

Although MERS-CoV has been shown to replicate in 
various animal cell lines including bat cells _ 8 , the broad 
tissue tropism was mainly demonstrated using type strain 
EMC/2012, and no comparison was made between dif¬ 
ferent MERS-CoV strains. Moreover, cells from the bat 
hosts of MERS-CoV-related viruses, such as T pachypus 
and P. abramus which harbor Ty-BatCoV HKU4 and Pi- 
BatCoV HKU5, respectively, were not included in pre¬ 
vious studies, which may be due to the geographical 
limitation of these bat species. To better understand the 
replicative ability of MERS-CoV in bat cells, which may 
provide clues on the origin of MERS-CoV, we developed 
diverse primary bat cell lines from different bat species, 
including Rhinolophus sinicus (the host of SARSr-Bat- 
CoV) and T. pachypus (the host of Ty-BatCoV HKU4), 


and tested their susceptibilities to infection by different 
strains of MERS-CoV, SARS-CoV, and HCoV-229E. The 
DPP4 mRNA sequences of six bat species and their 
expression in bat cells were determined to correlate with 
viral replication results. Our findings showed differential 
cell tropism between different strains of MERS-CoV, 
SARS-CoV, and HCoV-229E, which offers insights to the 
origin of MERS-CoV. 

Results 

Five of 12 tested bat cell lines are susceptible to MERS-CoV 
EMC/2012 infection 

Since lineage C betacoronaviruses closely related to 
MERS-CoV were detected in bats, we developed 12 diverse 
primary bat cell lines from seven different bat species, 
including Hipposideros pomona } Miniopterus pusillus , 
Myotis ricketti , Pipistrellus abramus (the host of Pi-BatCoV 
HKU5), Rhinolophus sinicus (the host of SARSr-BatCoV 
and Rs-BatCoV HKU2), Tylonycteris pachypus (the host of 
Ty-BatCoV HKU4), Rousettus leschenaultii (the host of 
many viruses including Ro-BatCoV HKU9), which were 
subjected to infection with MERS-CoV at multiplicity of 
infection (MOI) of 1. Viral titers were determined by RT- 
qPCR on day 5 p.i. Five of the 12 cell lines (M. ricketti lung, 
P. abramus kidney, R. sinicus kidney and lung, and R. 
leschenaulti kidney cells) and Vero cells propagated MERS- 
CoV with at least one logi 0 increase in viral load. The 
highest increase in viral load was observed in R. sinicus 
kidney and lung cells, which was comparable with that 
observed in Vero cells (Fig. 1). Cytopathic effects (CPEs) 
were observed in infected M. ricketti lung and R. sinicus 
lung cells with rounding of cells (Fig. 2). The infectivities of 
the viruses from culture supernatants were confirmed by 
passage in Vero cells with CPE. H. pomona kidney, M. 
pusillus kidney, M. ricketti kidney, P. abramus lung, R. 
leschenaulti lung, and T. pachypus kidney and lung cells did 
not support MERS-CoV infection. 

Different MERS-CoV strains displayed different infectivities 
on bat cells 

MERS-Co Vs are currently classified into three major 
clades, clade A, B, and C, which were further divided into 
subclades A1-A2, B1-B6, Cl, and non-Cl _52 . To test if 
different MERS-CoV strains may show similar infectivities 
on bat cells, the five bat cells (M. ricketti lung, P. abramus 
kidney, R. sinicus kidney and lung, and R. leschenaulti 
kidney cells), which were susceptible to MERS-CoV EMC/ 
2012 (belonging to clade Al), were subjected to infection by 
three other MERS-CoV strains isolated from camels in 
Dubai, D998/15 (belonging to clade A2, but with spike gene 
belonging to clade C), and D1189.1/15 and D1271.1/15 
(both belonging to clade B3) \ Among the five challenged 
bat cells, M. ricketti lung supported infection by both 
D 1189.1/15 and D1271.1/15; while R. sinicus lung and R. 
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Fig. 1 The twelve bat cell lines and Vero cells were subject to infection by MERS-CoV in clade A and clade B. The 12 bat cell lines (PAK: 
Pipistrellus abramus kidney, PAL Pipistrellus abramus lung, RSK: Rhinolophus sinicus kidney, RSL: Rhinolophus sinicus lung, MRK: Myotis ricketti kidney, 
MRL: Myotis ricketti lung, TPK: Tylonycteris pachypus kidney, TPL: Tylonycteris pachypus lung, HPK: Hipposideros pomona kidney, RLK: Rousettus 
leschenaultii kidney, RLL: Rousettus leschenaultii lung, MPK: Miniopterus pusillus kidney) and Vero cells were subject to infection by MERS-CoV EMC/ 
2012 (belonging to clade Al) with MOI of 1 (a). Culture supernatants were harvested at day 0 and 5 postinfection. Viral titers were determined by real¬ 
time quantitative RT-PCR. Viral load was expressed as log 10 copies/mL. Error bars indicate the standard deviation of triplicate samples. The five bat cell 
lines susceptible to MERS-CoV EMC/2012 infection with > 1 log 10 increase in viral load at day 5 were marked with red triangles. They were subject to 
infection by three other MERS-CoV strains isolated from camels in Dubai, D998/15 (belonging to clade A2) (b), D1189.1/15 (c) and D1271.1/15 
(d) (belonging to clade B3) (b). Different MERS-CoV strains displayed different infectivities in these five bat cells. (*P< 0.05; **P<0.01; ***P< 0.001) 


leschenaulti kidney supported infection by D1189.1/15, with 
at least one logi 0 increase in viral load. None of the cells 
supported infection by D998/15 (Fig. 1). 

SARS-CoV can replicate in R. sinicus cells 

Since Chinese horseshoe bats are the major reservoir of 
SARS-related-CoVs, we also tested the 12 bat cell lines, 
including R. sinicus cells, for susceptibility to a clinical 
strain of SARS-CoV. SARS-CoV strain HKU-39849 can 
replicate in R. sinicus kidney but not lung cells, with at 
least one log 10 increase in viral load (Fig. 3a). The other 
bat cells did not support SARS-CoV infection. 

HCoV-229E cannot replicate in bat cells 

The 12 bat cells were also tested for susceptible to HCoV- 
229E infection. HCoV-229E strain ATCC VR-740, pre¬ 
viously isolated from a man with upper respiratory illness, 
cannot replicate in any of the tested bat cells (Fig. 3b). 

mRNA transcript sequence analysis and expression of 
DPP4 in bat cells 

Partial DPP4 mRNA transcript sequences (corre¬ 
sponding to nt 688-1040 of hDPP4 which includes 


residue 229-346, where the critical residues for MERS- 
CoV spike protein binding are present) were determined 
for six of the seven bat species from which the 12 tested 
bat cells were developed. The sequence of M. pusillus was 
not determined, as RT-PCR for the DPP4 mRNA from bat 
cells was negative. Phylogenetic analysis showed that the 
bat DPP4 mRNA sequences formed a distinct cluster 
being closely related to sequences from human and non¬ 
human primates; while the dromedary camel 
DPP4 sequence was most closely related to that of wild 
boar (Fig. 4a). Previous studies have identified 14 critical 
residues in hDPP4 for binding of MERS-CoV spike pro¬ 
tein 3,54 Upon multiple sequence alignment of the cor¬ 
responding regions that contain these critical residues, 
most of the residues are conserved in the six bat 
DPP4 sequences in this study (Fig. 3b). Notably, T. 
pachypus DPP4 (Tp-DPP4) contains an I I< substitution 
at position 295 compared with hDPP4. On the other hand, 
residue R336 of hDPP4 was only conserved in T. pachypus 
and H. pomona, both not susceptible to MERS-CoV 
infection. 

RT-qPCR of bat DPP4 mRNA in the corresponding bat 
cells was performed to determine the mRNA expression 
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Fig. 2 Cytopathic effects (CPE) were observed in infected M. ricketti lung, /?. sinicus lung and Vero cells on 5 days postinfection. CPE was 

compared between Myotis ricketti lung (immortalized) cells that were uninfected (control) (a), and infected with Dubai camel MERS strains D998/15 
(b), D1189.1/15 (c) & D1271.1/15 (d). CPE was compared between Rhinolophus sinicus lung (immortalized) cells that were uninfected (control) (e), and 
infected with Dubai camel MERS strains D998/15 (f), D1189.1/15 (g), and D1271.1/15 (h). CPE was compared between Vero cells that were uninfected 
(control) (i), and infected with Dubai camel MERS strains D998/15 (g), D1189.1/15 (k), and D1271.1/15 (I) 
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Fig. 3 The twelve bat cell lines and Vero/HFL cells were subject to infection by SARS-CoV and HCoV-229E. The 12 bat cell lines (PAK: 
Pipistrellus abramus kidney, PAL Pipistrellus abramus lung, RSK: Rhinolophus sinicus kidney, RSL: Rhinolophus sinicus lung, MRK: Myotis ricketti kidney, 
MRL: Myotis ricketti lung, TPK: Tylonycteris pachypus kidney, TPL: Tylonycteris pachypus lung, HPK: Hipposideros pomona kidney, RLK: Rousettus 
leschenaultii kidney, RLL: Rousettus leschenaultii lung, MPK: Miniopterus pusillus kidney) and Vero/HFL cells were subject to infection by SARS-CoV 
with MOI of 1 (a) and HCoV-229E with MOI of 0.01(b). Culture supernatants were harvested at day 0 and 5 postinfection. Viral titers were determined 
by real-time quantitative RT-PCR. Viral load was expressed as log 10 copies/mL. Error bars indicate the standard deviation of triplicate samples. Only 
RSK cells can support SARS-CoV infection with > 1 log 10 increase in viral load at day 5 (blue triangle) and none of the 12 bat cell lines support HCoV- 
229E infection. (TP <0.05; **P< 0.01; ***P< 0.001) 


levels. Results showed that DPP4 mRNA was expressed in 
all the five bat cells that were susceptible to MERS-CoV 
infection, while it was also expressed in H . pomona kidney 
and R. leschenaulti lung cells which were not susceptible 
to MERS-CoV infection. The mRNA expression level in P. 


abramus kidney (susceptible to MERS-CoV) was sig¬ 
nificantly higher than that in its lung (non-susceptible) 
cells (P = 0.0185 by Students t test). Similarly, the mRNA 
expression level in R. leschenaulti kidney (susceptible to 
MERS-CoV) was significantly higher than that in its lung 
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Fig. 4 Phylogenetic analyses of partial DPP4 mRNA sequences and comparison of critical amino acid residues of human, camels, bats and 
other animals, (a) The trees were constructed by Neighbor-Joining method using substitution models and bootstrap values calculated from 
1000 trees. Only bootstrap values >70% are shown. One hundred and twelve aa positions were included in the analyses. The scale bars represent 
20 substitutions per site. Bat DPP4s that are sequenced in this study are labeled with black circles. Comparison of critical amino acid residues in DPP4 
from different animal host for receptor binding in the region of residues 229-346 with respect to human DPP4 (b) 
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Fig. 5 DPP4 expression analysis of bat cell lines, (a) mRNA levels of DPP4 were measured in various bat cells extracts by RT-qPCR and plotted 
relative to Vero cells, normalized by (3-actin mRNA levels. The mRNA levels of DPP4 were compared between susceptible and non-susceptible bat cell 
lines (b). Statistical significance was assessed by Student's f test P< 0.05 


(non-susceptible) cells (P = 0.0009 by Students t test). 
Compared with bat cells that are non-susceptible, bat cells 
that are susceptible to MERS-CoV showed a significantly 
higher mRNA expression level (P = 0.0174 by Students t 
test) (Fig. 5). 

Overexpression of Tp-DPP4 does not confer MERS-CoV 
EMC/2012 susceptibility to T. pachypus cells 

Since DPP4 was not expressed in T. pachypus cells 
while this bat species hosts Ty-BatCoV HKU4, a close 
relative of MERS-CoV which can utilize hDPP4 for cell 
entry 41,42 , we overexpressed Tp-DPP4 in T. pachypus lung 
and kidney cells for infection by MERS-CoV. While a 
slight increase (less than one log 10 ) in viral replication was 
noted in both mock and Tp-DPP4-overexpressed cells on 
day 5 p.i., no significant difference was noted between 
mock and Tp-DPP4-overexpressed cells (Fig. 6), sug¬ 
gesting that overexpression of Tp-DPP4 does not confer 
MERS-CoV susceptibility to T. pachypus cells. 

Discussion 

The present study provides further support to the bat 
origin of MERS-CoV. Previous studies have shown that 
MERS-CoV can infect cell lines from different animals, 
including cells from various bat species . However, bat 
cells from host species of MERS-CoV-related viruses, 
such as T. pachypus and P. abramus which harbor Ty- 
BatCoV HKU4 and Pi-BatCoV HKU5, respectively, were 
not included in these studies. In this study, MERS-CoV 
EMC/2012 was able to replicate in cells from four dif¬ 
ferent bat genera/species belonging to three different bat 
families including Pteropodidae , Rhinolophidae, and Ves- 
pertilionidae. None of these bat cells were previously 
reported to support MERS-CoV replication. While lineage 
C betacoronaviruses have not been detected in bats out¬ 
side the Vespertilionidae family, such diverse cellular 
tropism should prompt further studies to explore the host 
diversity of this group of CoVs and hence the possible 
evolutionary origin of MERS-CoV. In particular, the 


ability of MERS-CoV to replicate in cells from P. abra¬ 
mus , the host of Pi-BatCoV HKU5, may suggest that the 
ancestor of MERS-CoV could have originated from 
Pipistrellus- related bats. Interestingly, MERS-CoV 
showed the highest replicative ability in cells from R. 
sinicus . This suggests that MERS-CoV may potentially 
infect this bat species which is also the natural reservoir of 
SARSr-CoVs and animal origin of SARS-CoV 21 . Never¬ 
theless, three other MERS-CoV strains belonging to either 
clade A2 or B3 showed differential replicative abilities in 
the five bat cell lines that were susceptible to EMC/ 
2012 strain (clade Al), with strain D998/15 (clade A2) 
unable to replicate in all five tested cell lines. Interestingly, 
it has been shown in a recent study that some clade C 
MERS-CoV strains from Nigeria and Burkina Faso repli¬ 
cated at a lower titre in Calu-3 cells 2 . This is in line with 
the present study that our clade A2 strain D998/15 (with 
its spike gene belonging to clade C), conferred less broad 
cellular tropism. Our findings suggest that different clades 
or strains of MERS-CoV may possess different cellular 
tropism and host range, while all MERS-CoVs can repli¬ 
cate well in the immune-deficient Vero cells. 

The broader cellular tropism of MERS-CoV than SARS- 
CoV and HCoV-229E in bat cells may reflect the host 
diversity of lineage C betacoronaviruses in bats. Bats are 
now known to be the recent origin of at least two human 
Co Vs, including SARS-CoV and HCoV-229E. SARS-CoV 
is likely a recombinant virus arising from horseshoe bats 
before it jumped to civet as the intermediate host and 
then human 6,24 . SARSr-CoVs in bats were shown to 
utilize the SARS-CoV receptor, human angiotensin¬ 
converting enzyme 2 (hACE2), for cell entry 55 . Similar 
to SARS-CoV, HCoV-229E is likely to have originated 
from bats. CoVs closely related to HCoV-229E have been 
detected in bats of the genus Hipposideros in Africa 6 . 
More recently, CoVs even closer to HCoV-229E were 
identified in dromedary camels in the Middle East, which 
were able to utilize the receptor of HCoV-229E, human 
aminopeptidase N (hAPN) for cell entry . This suggests 























Lau et al. Emerging Microbes & Infections (2018)7:209 


Page 7 of 11 



that dromedary camels may have served as intermediate 
hosts for bat-to-human transmission of HCoV-229E, 
which may be analogous to the evolution of MERS-CoV. 
Recently, a CoV, closely related to HCoV-NL63 in most 
genome regions except the spike protein, was detected 
from a bat of the genus Triaenops in Kenya 58 . However, 
further studies are required to identify even closer rela¬ 
tives of HCoV-NL63 to ascertain its possible bat origin. In 
contrast to SARSr-CoVs and HCoV-229E which were 
mainly found in horseshoe and roundleaf bats, respec¬ 
tively, a more diverse host range was observed in lineage 
C betacoronaviruses. Yet, the different bat species har¬ 
boring lineage C betacoronaviruses all belonged to the 
family Vespertilionidae . Specifically, Ty-BatCoV HKU4, 
Pi-BatCoV HKU5, Hp-BatCoV HKU25, and BtVs-Beta- 
CoV/SC2013 were detected in bats belonging to the 
genera, Tylonycteris , Pipistrellus , Hypsugo, and Verspeti- 
lio , respectively, in China, whereas NeoCoV and BatCoV 
PREDICT/PDF-2180 were detected in bats belonging to 

O/T on 

Neoromicia and Pipistrellus in Africa ‘ ’ . In this study, 
SARS-CoV could only replicate in R. sinicus kidney cells 
but not other bat cells, which reflects its evolutionary 
origin from horseshoe bats. On the other hand, HCoV- 
229E was unable to replicate in any tested bat cells 
including cells from Hipposideros . This may suggest a 
relatively narrow host range of SARS-CoV and HCoV- 
229E in bats compared with lineage C betacoronaviruses 
including MERS-CoV. 

The close phylogenetic relationship between bat and 
primate DPP4 sequences may reflect the replicative 
ability of MERS-CoV in bat cells. Moreover, cells that 
supported MERS-CoV replication showed significantly 
higher DPP4 mRNA expression than those that were 


non-susceptible to MERS-CoV, suggesting that cellular 
DPP4 expression is critical for viral infection. Since Ty- 
BatCoV HKU4 from T. pachypus was previously shown 
to be able to utilize hDPP4 for cell entry 41,42 , we 
expected MERS-CoV to be able to utilize Tp-DPP4 for 
receptor binding and infect T. pachypus cells. While the 
inability of MERS-CoV to replicate in T. pachypus cells 
may be partly explained by the lack of DPP4 expression, 
overexpression of Tp-DPP4 was unable to confer viral 
susceptibility. Since mRNA expression may not always 
correlate with cell surface expression of the receptor, 
further studies are required to study the receptor¬ 
binding interphase between MERS-CoV and Tp-DPP4 
and whether other cellular factors may play a role in 
determining viral replicative abilities in T. pachypus 
cells, which may offer further insights into the evolu¬ 
tionary origin and mechanisms of interspecies trans¬ 
mission of MERS-CoV. Interestingly, one of the critical 
residues for binding to MERS-CoV spike protein, R336, 
was only conserved in T. pachypus and H. pomona DPP4 
(cells from these two species did not support MERS- 
CoV replication), but not in the other four sequenced 
bat DPP4 (cells from these four species supported 
MERS-CoV replication). This suggests that this residue 
may not be critical for receptor binding. Binding and 
mutagenesis studies may help better understand the role 
of receptor-binding-interphase during viral evolution 
and interspecies jumping. 

Materials and methods 

Ethics statement 

The collection of bat samples for developing bat cell 
lines was approved by the Committee on the Use of Live 
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Table 1 Cell lines used in the present study 


Organism 


Site of origin 

Source 

Bat 

Hipposideros pomona (Pomona roundleaf bat) 

Kidney 


n-house development 


Miniopterus pusillus (Lesser bent-winged bat) 

Kidney 


n-house development 


Myotis ricketti (Rickett's big-footed bat) 

Kidney 


n-house development 


Myotis ricketti (Rickett's big-footed bat) 

Lung 


n-house development 


Pipistrellus abramus (Japanese pipistrelle) 

Kidney 


n-house development 


Pipistrellus abromus (Japanese pipistrelle) 

Lung 


n-house development 


Rhinolophus sinicus (Chinese horseshoe bat) 

Kidney 


n-house development 


Rhinolophus sinicus (Chinese horseshoe bat) 

Lung 


n-house development 


Tylonycteris pochypus (Lesser bamboo bat) 

Kidney 


n-house development 


Tylonycteris pochypus (Lesser bamboo bat) 

Lung 


n-house development 


Rousettus leschenoultii (Leschenault's rousette) 

Kidney 


n-house development 


Rousettus leschenoultii (Leschenault's rousette) 

Lung 


n-house development 

Human 

Homo sapiens 

Embryonic Lung 

HFL (In-house development) 

Monkey 

Cercopithecus oethiops (African green monkey) 

Kidney 

Vero (ATCC CCL-81) 


Animals for Teaching and Research, the University of 
Hong Kong, Hong Kong, China. 

Cell lines 

The bat cells, Vero (African green monkey kidney) 
cells, and HFL (human embryonic lung fibroblast) cells 
used in this study are described in Table 1. For devel¬ 
opment of primary bat cell lines, bats were captured in 
Hong Kong and euthanized before dissection for 
recovery of cells from organs aseptically. Briefly, the 
kidney and lung tissue were rinsed with cold PBS and 
cut into small pieces. Cold 0.25% trypsin-EDTA was 
added to the tissues and incubated at 4 °C overnight. 
Tissues were then incubated at 37 °C on shaking plat¬ 
form for 30 min. Supernatants were filtered through cell 
strainers to remove large pieces of tissues 9 . Bat cells 
were harvested by spinning down the supernatant at 
1200 rpm for 8 min and were grown in DMEM/ 
F12 supplemented with 15% FBS. Vero and HFL cells 
were grown in MEM supplemented with 10% FBS. All 
cells were incubated at 37 °C with 5% C0 2 . 

Virus isolates 

MERS-CoV strain EMC/2012 was provided by Fouchier 
and Zaki et al. 11 . MERS-CoV strains D998/15, D1189.1/ 
15, and D1271.1/15 were isolated from dromedary camels 
in Dubai’ 1 . SARS-CoV strain HKU-39849 was isolated 
from the brother-in-law of the index patient in Hong 
Kong during the SARS epidemic ’. The HCoV-229E strain 
ATCC VR-740 was used. MERS-CoV isolates and SARS- 
CoV were propagated in Vero cells at MOI of 0.01 in 


MEM supplemented with 1% FBS. HCoV-229E was pro¬ 
pagated in HFL. 

Infection of bat cell lines 

Viral titers were determined as median tissue culture 
infective dose (TCID 50 ) per ml in confluent cells in 96- 
well microtiter plates. Cells were seeded onto 24-well 
tissue culture plates, at 2 x 10 5 cells per well with the 
respective medium and incubated at 37 °C and 5% C0 2 for 
24 h prior to experiment. Cells were washed once with 
PBS and inoculated with 1 MOI of MERS-CoV or SARS- 
CoV, or 0.01 MOI of 229E for 1 h. After 1 h of viral 
adsorption, the supernatant was removed and cells were 
washed twice with PBS. The cells were maintained in 
MEM supplemented with 1% FBS for Vero and HFL cells 
and DMEM/F12 supplemented with 1% FBS for bat cells, 
before further incubation for 5 days. 

Viral replication studies 

To study viral replication efficiency, progeny viruses 
from cell culture supernatants collected at 5 days post¬ 
infection (p.i.) were subjected to reverse transcription- 
quantitative PCR (RT-qPCR) according to our previous 
protocols °. Briefly, total RNA extracted from cell 
culture supernatants with QIAsymphony DSP Virus/ 
Pathogen Mini Kit (Qiagen, Hilden, Germany) was 
reverse transcribed and amplified with MERS-CoV 
primers (forward primer 5'-CAAAACCTTCCCTAAGA 
AGGAAAAG-3' reverse primer 5'-GCTCCTTTGGAG 
GTTCAGACAT-3'), SARS-CoV primers (forward 5'- AC 
CAGAATGGAGGACGCAATG-3' reverse 5'-GCTGTG 
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Table 2 Primers used for RT-qPCR in this study 


Target 

Primers 




Forward 

Reverse 

Probe 

MERS- 

CoV N 

gene 

5'-CAAAACCTTCCCTAAGAAGGAAAAG 

-3' 

5'-GCTCCTTT GG AGGTT CAG ACAT-3' 

(5'-(FAM)ACAAAAGGCACCAAAAGAAGAATCAAC 

AGACC(BHQ1)-3' 

SARS- 

CoV N 

gene 

5'- ACCAGAATGGAGGACGCAATG-3' 

5'-GCTGTGAACCAAGACGCAGTATTAT-3' 

5'-(FAM)ACCCCAAGGTTTACCC(NFQ)-3' 

HCoV- 

229E N 

gene 

5'-CAGTCAAATGGGCTGATGCA-3' 

5'-AAAGGGCTATAAAGAGAATAAGGTATTCT-3' 

5'-(FAM)CCCTGACGACCACGTTGTGGTrCA(BHQ1)- 

3' 

DPP4 

5'-TGATCTTGCCTCCTCATmGATAA-3' 

5'-GTAACCACTOCTCTGCCATCAA-3' or 5'- 

GTAACCAC^CCTCTGCCGTCAA-3' (for Rousettus 

cell lines) 

5'-(FAM)CCACMTTCAMACTCARYTGGGCTACTTAC 

C(BHQ1)-3' 

(3-actin 

5'-CTCTTCCAGCCCTCCTTCCT-3' 

S'-TTCATCGTGCTGGGAGCC-S' or 5'- 

TTCATTGTGCTGGGAGCC-3' (for Rousettus cell 

lines) or 5'-TPCATGGTGCTGGGGGCC-3' (for 
Rhinolophus cell lines) 

5'-(F AM)CAT G AAGT GYG ACGTBG AC AT CCG 

(BHQD-3' 


AACCAAGACGCAGTATTAT-3') and HCoV-229E pri¬ 
mers (forward 5 , -CAGTCAAATGGGCTGATGCA-3 / 
reverse 5'-AAAGGGCTATAAAGAGAATAAGGTAT 
TCT-3') using real-time one-step quantitative RT-PCR 
assay as described previously with modifications 5,6 °. 
Probes for MERS-CoV [5' - (F AM) AC AAAAGGC AC 
CAAAAGAAGAATCAACAGACC(BHQl)-3'], SARS- 
CoV [5 , -(EAM)ACCCCAAGGTTTACCC(NFQ)-3 / ] and 
HCoV-229E [5' - (F AM) CCCT G ACG ACC ACGTTGTG 
GTTCA(BHQl)-3'] were used (Table 2). Reactions were 
first incubated at 50 °C for 30 min, followed by 95 °C for 2 
min, and were then thermal cycled for 50 cycles (95 °C for 
15 s, 55 °C for 30 s). A series of 10 logi 0 dilutions 
equivalent to 1 x 10 1 to 1 x 10 10 copies per reaction 
mixture were prepared to generate calibration curves and 
were run in parallel with the test samples. All experiments 
were performed in triplicate. 

Amplification and sequencing of partial bat DPP4 mRNA 
transcripts 

Total RNA was extracted from bat cell lysates using 
RNeasy Mini Spin Column (QIAgen). cDNA was PCR 
amplified with primers, 5'-GTCACCAGAGGGTCAT 
AAA-3' and S'-CCACTTCCTCTGCCATCAAA-S'. The 
PCR mixture (25 pi) contained cDNA, PCR buffer, 200 
pM (each) dNTPs, and 1.0 U Iproof Polymerase (Bio-Rad, 
Hercules, CA, USA). The mixtures were amplified for 40 
cycles of 98 °C for 10 s, 55 °C for 30 s, and 72 °C for 72 s 
and a final extension at 72 °C for 10 min in an automated 


thermal cycler (Applied Biosystems, Foster City, CA, 
USA). RT-PCR products were gel purified using QIA- 
quick gel extraction kit (Qiagen), and sequenced with an 
ABI Prism 3700 DNA Analyzer (Applied Biosystems). The 
sequences obtained were compared with sequences of 
DPP4 genes in GenBank database. Phylogenetic tree 
construction was performed based on an amino acid 
alignment of partial DPP4 sequences (corresponding to 
residue 229-346 of hDPP4) using the Neighbor-Joining 
method with JTT model by MEGA 6.0, with bootstrap 
values calculated from 1000 trees. 

DPP4 expression analysis 

To study DPP4 expression profiles in different bat cell 
lines, cell lysates were collected for total RNA extraction 
using RNeasy Mini Spin Column (QIAgen). RNA was 
eluted in 50 pi of RNase-free water and was used as a 
template for one-step RT-qPCR with Superscript III 
platinum One-step qRT-PCR system (Invitrogen, Carls¬ 
bad, CA, USA). RT-qPCR assays were performed using 
conserved primers designed by multiple alignment of 
available bat DPP4 gene sequences (Table 2), using 
fi-actin for normalization. cDNA was amplified in a 
EightCycler 480 (Roche, Basel, Switzerland) with 25 pi 
reaction mixtures containing 2 x reaction mix, 5 pi RNA, 
25 pM ROX reference dye, 50 pM primers, and 10 pM 
probe at 50 °C for 30 min, then 95 °C for 2 min followed 
by 50 cycles of denaturation, annealing and extension. 
Experiments were performed in triplicates, and results 
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were expressed as the mean expression level of DPP4/ 
fi-actin . The relative expression between different bat cells 
was then calculated by the AACt method. 

tpDPP4 overexpression in T. pachypus cells 

T. pachypus DPP4 (tpDPP4) sequence was cloned into 
pLenti7.3/V5-TOPO vector (Invitrogen). The construct 
was transfected into 293FT cells together with ViraPower 
Packaging Mix (Invitrogen) using Lipofectamine 2000 
(Life Technologies, Carlsbad, CA, USA). Lentivirus was 
harvested from the supernatant and concentrated using 
PEG-it (System Biosciences, Palo Alto, USA). T. pachypus 
cells resistant to MERS-CoV were transduced using the 
concentrated lentivirus for tpDPP4 overexpression and 
were subsequently subjected to inoculation with 1 MOI of 
MERS'CoV. 

Nucleotide sequence accession numbers 

The nucleotide sequences of bat DPP4 obtained from 
this study have been deposited in the GenBank sequence 
database under accession numbers MH345671- 
MH345676. 

Acknowledgements 

We thank C.T. Shek, Agriculture, Fisheries, and Conservation Department, the 
Government of the HKSAR, for capture of bats for bat cell development and 
expert opinion. This work is partly supported by Theme-Based Research 
Scheme (T11-707/15-R) and Research Grant Council Grant, University Grant 
Council; Health and Medical Research Fund of the Food and Health Bureau of 
HKSAR; University Development Fund and Special Research Achievement 
Award, The University of Hong Kong; Consultancy Service for Enhancing 
Laboratory Surveillance of Emerging Infectious Disease for the HKSAR 
Department of Health; The Collaborative Innovation Center for Diagnosis and 
Treatment of Infectious Diseases, the Ministry of Education of China; National 
Science and Technology Major Project of China (grant number 
2012ZX10004213). Views expressed in this paper are those of the authors only, 
and may not represent the opinion of the AFCD or the Government of the 
HKSAR. 

Author details 

1 State Key Laboratory of Emerging Infectious Diseases, Li Ka Shing Faculty of 
Medicine, The University of Hong Kong, Hong Kong, China, department of 
Microbiology, Li Ka Shing Faculty of Medicine, The University of Hong Kong, 
Hong Kong, China. 3 Carol Yu Centre for Infection, Li Ka Shing Faculty of 
Medicine, The University of Hong Kong, Hong Kong, China. Collaborative 
Innovation Center for Diagnosis and Treatment of Infectious Diseases, Li Ka 
Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China. 
Central Veterinary Research Laboratory, Dubai, United Arab Emirates 

Author contributions 

S.K.P.L. and P.C.Y.W. conceived the study and designed the experiments. R.Y.Y.F., 

H.K.H.L., L.C., J.F., K.S.M.L., E.Y.M.W. and S.S.A. performed the experiments. S.K.P.L., 
R.F.F.Y., H.K.H.L., J.F., K.S.M.L., J.F.W.C., R.K.H.K., K.H.C., U.W., K.Y.Y. and P.C.Y.W. 
analyzed the data. S.K.P.L., R.F.F.Y., H.K.H.L., K.S.M.L., L.C. and P.C.Y.W. wrote the 
paper. 

Conflict of interest 

The authors declare that they have no conflict of interest. 


Publisher's note 

Springer Nature remains neutral with regard to jurisdictional claims in 
published maps and institutional affiliations. 


Received: 24 August 2018 Revised: 8 November 2018 Accepted: 11 
November 2018 

Published online: 10 December 2018 


References 

1. de Groot R. J. et al. Coronaviridae. In: Virus Taxonomy, Classification and 
Nomenclature of Viruses. Ninth Report of the International Committee on 
Taxonomy of Viruses, International Union of Microbiological Societies, Virology 
Division, 806-828 (2011). 

2. Woo, P. C. et al. Discovery of seven novel Mammalian and avian coronaviruses 
in the genus deltacoronavirus supports bat coronaviruses as the gene source 
of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene 
source of gammacoronavirus and deltacoronavirus. J. Virol. 86 , 3995-4008 
( 2012 ). 

3. Woo, P. C. et al. Comparative analysis of twelve genomes of three novel group 
2c and group 2d coronaviruses reveals unique group and subgroup features. 

J. Virol. 81 , 1574-1585 (2007). 

4. Gorbalenya, A. E., Snijder, E. J. & Spaan, W. J. M. Severe acute respiratory 
syndrome coronavirus phylogeny: toward consensus. J. Virol. 78 , 7863-7866 
(2004). 

5. Drosten, C. et al. Identification of a novel coronavirus in patients with severe 
acute respiratory syndrome. New Engl. J. Med. 348 , 1967-1976 (2003). 

6. Fouchier, R. A. M. et al. A previously undescribed coronavirus associated with 
respiratory disease in humans. Proc. Natl Acod. Sci. 101 , 6212-6216 (2004). 

7. Ksiazek, T. G. et al. A novel coronavirus associated with severe acute respiratory 
syndrome. New Engl. J. Med. 348 , 1953-1966 (2003). 

8. Peiris, J. S. M. et al. Coronavirus as a possible cause of severe acute respiratory 
syndrome. Lancet 361 , 1319-1325 (2003). 

9. van der Hoek, L. et al. Identification of a new human coronavirus. Nat. Med. 10 , 
368-373 (2004). 

10. Woo, P. C. et al. Characterization and complete genome sequence of a novel 
coronavirus, coronavirus HKU1, from patients with pneumonia. J. Virol. 79 , 
884-895 (2005). 

11. Zaki, A. M., van Boheemen, S., Bestebroer, T. M., Osterhaus, A. D. & Fouchier, R. 
A. Isolation of a novel coronavirus from a man with pneumonia in Saudi 
Arabia. New Engl. J. Med. 367 , 1814-1820 (2012). 

12. de Groot, R. J. et al. Middle East respiratory syndrome coronavirus 
(MERS-CoV): announcement of the Coronavirus Study Group. J. Virol. 87 , 
7790-7792 (2013). 

13. Woo, P. C. Y. et al. Comparative analysis of 22 coronavirus HKU1 genomes 
reveals a novel genotype and evidence of natural recombination in cor¬ 
onavirus HKU1. J. Virol. 80 , 7136-7145 (2006). 

14. Kottier, S. A., Cavanagh, D. & Britton, P. Experimental evidence of recombi¬ 
nation in coronavirus infectious bronchitis virus. Virology 213 , 569-580 (1995). 

15. Lau, S. K. P. et al. Isolation and characterization of a novel betacoronavirus 
subgroup a coronavirus, rabbit coronavirus HKU14, from domestic rabbits. J. 
Virol. 86 , 5481-5496 (2012). 

16. Lau, S. K. P. et al. Ecoepidemiology and complete genome comparison of 
different strains of severe acute respiratory syndrome-related rhinolophus 
bat coronavirus in china reveal bats as a reservoir for acute, self- 
limiting infection that allows recombination events. J. Virol. 84 , 2808-2819 
( 2010 ). 

17. Lau, S. K. P. et al. Molecular epidemiology of human coronavirus OC43 
reveals evolution of different genotypes over time and recent emergence of 
a novel genotype due to natural recombination. J. Virol. 85 , 11325-11337 
( 2011 ). 

18. Lau, S. K. P. et al. Coexistence of different genotypes in the same bat and 
serological characterization of rousettus bat coronavirus HKU9 belonging to a 
novel betacoronavirus subgroup. J. Virol. 84 , 11385-11394 (2010). 

19. Lau, 5. K. P. et al. Discovery of a novel coronavirus, China Rattus Coronavirus 
HKU24, from Norway rats supports the murine origin of betacoronavirus 1 and 
has implications for the ancestor of betacoronavirus lineage A. J. Virol. 89 , 
3076-3092 (2014). 

20. Guan, Y. Isolation and characterization of viruses related to the SAR5 
coronavirus from animals in southern China. Science 302 , 276-278 
(2003). 

21. Lau, S. K. et al. Severe acute respiratory syndrome coronavirus-like virus in 
Chinese horseshoe bats. Proc. Natl Acad. Sci. USA 102 , 14040-14045 (2005). 

22. Li, W. Bats are natural Reservoirs of SARS-like coronaviruses. Science 310 , 
676-679 (2005). 



Lau et al. Emerging Microbes & Infections (2018)7:209 


Page 11 of 11 


23. Lau, 5. K. et al. Severe acute respiratory syndrome (SARS) coronavirus ORF8 
protein is acquired from SARS-related coronavirus from greater horseshoe 
bats through recombination. J. Virol. 89 , 10532-10547 (2015). 

24. Song, H. D. et al. Cross-host evolution of severe acute respiratory syndrome 
coronavirus in palm civet and human. Proc. Natl Acad. Sci. 102 , 2430-2435 
(2005). 

25. Chinese, S. M. E. C. Molecular evolution of the SARS coronavirus during the 
course of the SARS epidemic in China. Science 303 , 1666-1669 (2004). 

26. Lau, S. K. P. et al. Recent transmission of a novel alphacoronavirus, bat cor¬ 
onavirus HKU10, from Leschenault's Rousettes to Pomona Leaf-Nosed Bats: 
first evidence of interspecies transmission of coronavirus between bats of 
different suborders. J. Virol. 86 , 11906-11918 (2012). 

27. Woo, P. C. et al. Discovery of a novel bottlenose dolphin coronavirus reveals a 
distinct species of marine mammal coronavirus in Gammacoronavirus. J. Virol. 
88 , 1318-1331 (2014). 

28. Zhang, J. et al. Genomic characterization of equine coronavirus. Virology 369 , 
92-104 (2007). 

29. Elasoksuz, M. et al. Biologic, antigenic, and full-length genomic characterization 
of a bovine-like coronavirus isolated from a giraffe. J. Virol. 81 , 4981-4990 
(2007). 

30. Jin, L. et al. Analysis of the genome sequence of an alpaca coronavirus. 
Virology 365 , 198-203 (2007). 

31. Reusken, C. B. et al. Middle East respiratory syndrome coronavirus neutralising 
serum antibodies in dromedary camels: a comparative serological study. 
Lancet Infect. Dis. 13, 859-866 (2013). 

32. Chan, J. F. W. et al. Middle East respiratory syndrome coronavirus: another 
zoonotic betacoronavirus causing SARS-like disease. Clin. Microbiol. Rev. 28 , 
465-522 (2015). 

33. Woo, P. C. Y. et al. Molecular diversity of coronaviruses in bats. Virology 351 , 
180-187 (2006). 

34. Lau, 5. K. et al. Genetic characterization of betacoronavirus Lineage C viruses in 
bats reveals marked sequence divergence in the spike protein of pipistrellus 
bat coronavirus HKU5 in Japanese pipistrelle: implications for the origin of the 
novel Middle East respiratory syndrome coronavirus. J. Virol. 87 , 8638-8650 
(2013). 

35. Yang, L. et al. MERS-related betacoronavirus in Vespertilio superans bats, China. 
Emerg. Infect. Dis. 20 , 1260-1262 (2014). 

36. Corman, V. M. et al. Rooting the phylogenetic tree of middle East respiratory 
syndrome coronavirus by characterization of a conspecific virus from an 
African bat. J. Virol. 88 , 11297-11303 (2014). 

37. Anthony, S. J. et al. Further evidence for bats as the evolutionary source of 
Middle East respiratory syndrome coronavirus. mBio 8 (2017). 

38. Geldenhuys, M. et al. A metagenomic viral discovery approach identifies 
potential zoonotic and novel mammalian viruses in Neoromicia bats within 
South Africa. PloS ONE 13, eOl 94527 (2018). 

39. Lau, S. K. P. et al. Receptor usage of a novel bat lineage C betacoronavirus 
reveals evolution of MERS-related coronavirus spike proteins for human DPP4 
binding. J. Infect. Dis. 218,197-207 (2018). 

40. Corman, V. M. et al. Characterization of a novel betacoronavirus related to 
middle East respiratory syndrome coronavirus in European hedgehogs. J. Virol. 
88 , 717-724 (2014). 


41. Yang, Y. et al. Receptor usage and cell entry of bat coronavirus EIKU4 provide 
insight into bat-to-human transmission of MERS coronavirus. Proc. Natl Acad. 
Sci. USA 111 , 12516-12521 (2014). 

42. Wang, Q. et al. Bat origins of MERS-CoV supported by bat coronavirus EIKU4 
usage of human receptor CD26. Cell host & Microbe 16, 328-337 (2014). 

43. Cai, Y. et al. CD26/DPP4 cell-surface expression in bat cells correlates with 
bat cell susceptibility to Middle East respiratory syndrome coronavirus 
(MERS-CoV) infection and evolution of persistent infection. PloS ONE 9 , 
el 12060 (2014). 

44. Munster, V. J. et al. Replication and shedding of MERS-CoV in Jamaican fruit 
bats (Artibeus jamaicensis). Sci. Rep. 6 , 21878 (2016). 

45. Chan, J. F. et al. Differential cell line susceptibility to the emerging novel 
human betacoronavirus 2c EMC/2012: implications for disease pathogenesis 
and clinical manifestation. J. Infect. Dis. 207 , 1743-1752 (2013). 

46. Eckerle, I. et al. Replicative capacity of MERS coronavirus in livestock cell lines. 
Emerg. Infect. Dis. 20 , 276-279 (2014). 

47. Vergara-Alert, J. et al. Livestock susceptibility to infection with Middle East 
respiratory syndrome coronavirus. Emerg. Infect. Dis. 23 , 232-240 (2017). 

48. Banerjee, A. et al. Generation and characterization of eptesicus fuscus (Big 
brown bat) kidney cell lines immortalized using the myotis polyomavirus large 
T-antigen. J. Virol. Methods 237 , 166-173 (2016). 

49. Lau, S. K. P., Wong, A. C. P., Lau, T. C. K. & Woo, P. C. Y. Molecular evolution of 
MERS coronavirus: dromedaries as a recent intermediate host or long-time 
animal reservoir? Int. J. Mol. Sci. 18 (2017). 

50. Gotten, M. et al. Transmission and evolution of the Middle East respiratory 
syndrome coronavirus in Saudi Arabia: a descriptive genomic study. Lancet 
382 , 1993-2002 (2013). 

51. Lau, S. K. et al. Polyphyletic origin of MERS coronaviruses and isolation of a 
novel clade A strain from dromedary camels in the United Arab Emirates. 
Emerg. Microbes Infect. 5 , el 28 (2016). 

52. Chu, D. K. W. et al. MERS coronaviruses from camels in Africa exhibit region- 
dependent genetic diversity. Proc. Natl Acad. Sci. USA 115 , 3144-3149 (2018). 

53. Lu, G. et al. Molecular basis of binding between novel human coronavirus 
MERS-CoV and its receptor CD26. Nature 500 , 227-231 (2013). 

54. Wang, N. et al. Structure of MERS-CoV spike receptor-binding domain com- 
plexed with human receptor DPP4. Cell Res. 23 , 986-993 (2013). 

55. Li, W. et al. Angiotensin-converting enzyme 2 is a functional receptor for the 
SARS coronavirus. Nature 426 , 450-454 (2003). 

56. Corman, V. M. et al. Evidence for an ancestral association of human cor¬ 
onavirus 229E with bats. J. Virol. 89 , 11858-11870 (2015). 

57. Corman, V. M. et al. Link of a ubiquitous human coronavirus to dromedary 
camels. Proc. Natl Acad. Sci. USA 113 , 9864-9869 (2016). 

58. Tao, Y. et al. Surveillance of bat coronaviruses in Kenya identifies relatives of 
human coronaviruses NL63 and 229E and their recombination history. J. Virol. 
91 (2017). 

59. Crameri, G. et al. Establishment, immortalisation and characterisation of pter- 
opid bat cell lines. PloS ONE 4 , e8266 (2009). 

60. Lau, S. K. et al. Delayed induction of proinflammatory cytokines and sup¬ 
pression of innate antiviral response by the novel Middle East respiratory 
syndrome coronavirus: implications for pathogenesis and treatment. J. Gen. 
Virol. 94 , (Pt 12) 2679-90 (2013). 



